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THE ESTIMATION OF THE DISPERSION OF 
WINDBORNE MATERIAL 
By F. PASQUILL, D.Sc. 


Introduction.— [he theoretical estimation of the concentrations arising trom 
sources of gaseous or finely divided particulate material has for long been based 
on treatments of atmospheric diffusion developed by Sir Graham Sutton, and 
subsequently expressed in the well known formulae for surface and elevated 
sources (see Micrometeorology, p. 2g0'). These formulae are reliable for specifying 
the average distribution, over a few hundred metres downwind of a source 
operating for a few minutes on level unobstructed terrain, with a steady wind 


direction and neutral conditions of atmospheric stability. Extension to other 
circumstances has depended on empirical and often speculative adjustments of 
the diffusion parameters. 


During the last few years, further investigations at the Chemical Defence 
Experimental Establishment, Porton, have shown that a fairly rational allow- 
ance can now be made for the effects of much of the wide variation in atmos- 
pheric turbulence which occurs in reality. This progress includes some extension 
to longer distances of travel. Although many aspects of the problem require 
further attention, these recent developments, supported also by experimental 
studies in the U.S.A., form a basis for a tentative system of estimating diffusion 
in a wide range of meteorological conditions and over distances up to about 
100 kilometres. 


The purpose of this article is to review the recent background of theoretical 
and experimental results, and to give details of the proposed system of calcu- 
lating the distribution of concentration downwind of a source. These details 
are set out in two appendixes, the first giving complete instructions for carrying 
out the calculations, the second presenting an example. 


The crosswind spread of a cloud from a continuous point source. 
In the earlier treatments of diffusion from a continuous point source the main 
feature was the adoption of a conjectured form for the Lagrangian correlation 
coefficient in G. I. Taylor’s (1921) statistical analysis of the dispersion of 
particles (see Micrometeorology, p. 284"). Justification of the particular form was 
provided by using it to deduce a correct expression for the aerodynamic drag 





of a surface. In effect this meant that a correct form for the eddy diffusivity 
determining the vertical transport of momentum had been specified. The 
remaining important steps were the assumptions of identity in the diffusivities 
for momentum and airborne material, and of analogy in the laws governing 
vertical and lateral diffusion from a source in the lower atmosphere. 


Taylor’s analysis assumes that the field of turbulence affecting the particles 
is homogeneous and steady, that is, that its statistical properties do not depend 
on position or time. It would be difficult, if not impossible, to find any example 
of atmospheric turbulence in which these conditions are strictly satisfied, and 
indeed it is known that near the ground the structure of turbulence changes 
systematically with height. However, subject to certain restrictions, it is not 
unreasonable to assume the existence of quasi-homogeneous, quasi-steady pro- 
perties in the horizontal plane, and even in the vertical plane well away from 
the ground. The obvious restrictions are that the structure of the flow should 
not be systematically patterned by dynamical or thermal! influences associated 
with surface topography, and that its properties should not change radically 
over the period of interest (hence, if there is a large diurnal! variation, for 
example, attention should be confined to small fractions of a day at a time). 
For such quasi-homogeneous, quasi-steady conditions a simple and direct 
adaptation of Taylor’s treatment has been developed by Hay and Pasquill?. 


The essential step in the new treatment is the adoption of a simple hypothesis 
regarding the Lagrangian variations of velocity (that is, those experienced by 
a single particle as it travels), which are difficult if not impossible to measure, 
and the variations which can be observed by using an instrument at a fixed 
position. In terms of the auto-correlation coefficient of eddy velocity the 
hypothesis is that this function has the same shape (with regard to time-lag) 
in the two cases, but that the Lagrangian coefficient takes 8 times as long as the 
“fixed-point” coefficient to decay to a given magnitude. It can be demonstrated 
analytically that identity in shape is not really critical, the important require- 
ment being that the ratio (8) of the corresponding integrals of the whole 
correlograms should be known. Even then it follows directly from Taylor’s 
original treatment that the magnitude of diffusion is insensitive to 8 at short 
range, and only depends on f* at long range. 


In the simple form used in practice the crosswind spread of particles from 
their mean position is given by 


Oy|x = [Gole, simp eee 


where oy, is the standard deviation of the crosswind displacements of the par- 
ticles at a distance x downwind, gy is the standard deviation of wind direction, 
and i is the mean wind speed. The subscripts are used to denote that this 
standard deviation is obtained by forming averages of the wind direction over 
moving intervals x/#8, and using the values so obtained over a duration rf, 
equal to the duration of release of the material, or to the duration of sampling 
(or exposure to) the cloud, whichever is the shorter. When x/i is greater than 
the duration of release (that is, when the plume is detached from the source), 
or than the duration of sampling of a continuous plume, the diffusion will be 
determined in a complex way by larger and larger eddies than those which 
contribute to the variation of wind direction over the time t, and equation (1) 
will then underestimate o,. 
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From measurements of the crosswind spread of particles at a distance of 100 
metres from a source of duration three minutes [tT=3 min in equation (1)], 
values of B were deduced which varied considerably but averaged about four. 
It was also shown that this value provided a satisfactory interpretation of 
diffusion data obtained earlier at Porton, and more recently in the U.S.A., 
for distances of travel up to about 1000 metres. Many more determinations of 
the effective value of 8 will be required to give a reasonably complete descrip- 
tion, but for many practical purposes it would appear that a useful range of 
conditions is adequately represented by the foregoing value. 

For ranges of travel much longer than 1000 metres the method is open to 
question on the grounds that as the plume spreads vertically its lateral spread 
will be affected by the systematic variation of wind direction with height above 
the ground. However, a limited examination (three cases) of the crosswind 
spread at a distance of about 75 kilometres showed that equation (1) did in fact 
give a reasonable approximation in terms of the surface wind. These cases were 
in daytime, with fairly vigorous mixing over a depth of about 1000 metres. The 
values of a,/x and o, were as follows, the latter being determined from a 
Baxendall wind-direction recorder at the site of release: 


Oy/x 0*073 0*077 0° 061 


[Gg], x 4ia 0° 084 0: 066 0-068 radians 


Vertical diffusion at short range from a source at ground level.—-A 
general treatment of vertical diffusion in the lowest layers of the atmosphere, 
based on laboratory laws relating wind profile and surface drag, has been given 
by Calder.* The full exploitation of the treatment requires careful specification 
of the vertical profile of mean wind velocity characteristic of the surface. For 
practical purposes, however, changes in the small-scale roughness (for example, 
in grass length) have only a minor effect on diffusion. For example (see Calder’s 
paper, p. 166%), a six-fold increase in the roughness parameter, corresponding 
to a doubling of the aerodynamic drag, produces only a 25 per cent increase in 
the Aeight of cloud at 100 metres from a source (defined here as the height at 
which the concentration in the cloud is one-tenth of the ground-level value). So 
although the treatment has only been verified for grassland, the implication is 
that increases in roughness of a moderate order (for example, due to crops, 
smal] bushes and hedges) are unlikely to make a vast difference to the vertical 
diffusion. This suggests that for open country, but away from larger disturbances 
such as those caused by woods, buildings or sharp changes in contour, an 
acceptable working approximation is probably best represented by the long- 
grass case (see Calder’s paper, p. 166%), for which the heights of cloud at 100 
and 1000 metres downwind are respectively 10 and 70 metres. 


It should be emphasized that these figures apply only to neutral conditions 
of stability, and cannot be directly extrapolated to longer ranges. The latter 
restriction arises from the assumption in the theoretical argument that the 
horizontal shearing stress is constant with height. In the present state of know- 
ledge this can be assumed with confidence only over the first few tens of metres. 

The method has been extended to non-neutral conditions by Deacon,‘ using 
an empirical power-law form of wind profile which reduces to the required 
logarithmic form in the special case of neutral conditions. There is evidence for 
reasonable agreement with observations of diffusion in unstable conditions, 
though the accuracy required in measuring the wind profile is even greater 





than in neutral conditions and there are other difficulties when the roughness 
elements are easily bent over or distorted by the wind. In stable conditions the 
method gives discrepant results, and Deacon® has recently argued that this may 
be because of the flow not being aerodynamically rough at the very low wind 
speeds which then occur very close to the ground. More recently, treatments 
by Monin® have attracted much attention from an analytical point of view, but 
have yet to be adapted for practical use. 


At the present stage the most useful guide to the estimation of vertical spread 
in non-neutral conditions is provided by the Prairie Grass measurements: in the 
U.S.A., which have been summarized in a convenient form by Cramer.’ These 
results show that cloud height at 100 metres downwind of the source ranged from 
about 4 metres in “‘extreme stability’’ to about 25 metres in “‘extreme instability”, 
with 74 metres in neutral conditions. As the site was relatively smooth (rough- 
ness parameter 2) < 1 cm) the latter figure compares satisfactorily with the value 
of 8-1 calculated by Calder (loc. cit.) for a z, of 0-5 centimetres. The magnitudes 
of vertical spread at greater distances (up to 800 metres) were inferred by 
Cramer from the measured variations (with distance) of peak concentration 
and cloud width. In unstable conditions they show an acceleration of the vertical 
spread with increasing distance. For example, in ‘“‘extreme instability” the 
inferred cloud height at 800 metres downwind is about 1000 metres. Because 
of the indirect derivation (including the assumption of Gaussian shape in the 
vertical distibution) it is difficult to judge the precision of this latter value. 
However, such a rate of spread merely requires the incidence of sustained 
upcurrents with an inclination of about 45°, and there is little doubt that such 
upcurrents do occur with light winds and well developed convection. 


The effect of an elevation of the source.—When material is released at 
an elevated position it may be expected that the effects of the variation of wind 
structure with height will initially be of secondary importance. As a first 
approximation, therefore, the vertical spread at short range may be derived on 
the assumption of quasi-homogeneous turbulence. In this case, on the same 
lines as those followed for lateral spread, the standard deviation of vertical 
spread (¢,) at distance x is given by 


ae/x = [O4]s, Tis er 


where o¢ is the standard deviation of the wind inclination (in radians). The 
latter is obtained from data averaged over periods 7/8 (7 =x/ia), and observed 
over a period t equal to the duration of release (or of sampling) of the material. 
The application of equation (2) has recently been tested against experimental 
data obtained in the U.S.A. 


Following custom, the plume from an elevated source is assumed to have 
Gaussian distributions of material both laterally and vertically, with standard 
deviations a, and o;. For simplicity, wind velocity is taken to be constant with 
height, and it is assumed that the effect of the ground can be represented by an 
image source (as in Micrometeorology, pp. 139 and 292"). With the usual co-ordinate 
system (x alongwind, y crosswind, z vertical) the continuity condition leads to 
the following expression for the distribution of concentration y (x, y, z) from a 
continuous source of strength Q at position (0, 0, H). 


__ 98 2 f —= fd 2 H 2) 
1 (89,2) = SAP es) [exp | + exp{ S20 | + 


220,0,0 20," 20," 
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With the simplifying assumption that o,/¢, is a constant independent of distance, 
differentiation of the above expression with respect to a, leads to the condition 
that y (x, 0, o) (that is, the concentration on the axis of the plume at ground 
level) is a maximum when 


0, V2 =H. 


+ 00 
Similarly, if the integral of the crosswind distribution, that is | 7 (x,y, 0)dy, 
oo 


is considered, the condition for this quantity to be a maximum is 
CG, = fi. oso oe 
Measurements have recently been made at the National Reactor Testing 
Establishment, Idaho Falls, U.S.A., of the ground-level distribution of a 
fluorescent tracer released at a height of about 50 metres. These provided 
observed estimates of the distance, d (max), at which the crosswind integrated 
concentration was a maximum. Records were also taken of the fluctuating 
inclination of the wind near the point of release, from which it was possible to 
evaluate o,, for specified values of t and x/a8 as required by equation (2). 
Substituting the condition (5) in equation (2) gives 


H = d (max) [0,],,, cos 


with s=d (max)/a8, from which calculated estimates of d (max) were obtained 
(assuming 8=4), the observed estimates being used in assigning the appropriate 
averaging time for analysing the wind trace. The observed and calculated 
estimates from thirteen tests are listed in Table |. 


TABLE I— OBSERVED AND CALCULATED ESTIMATES OF THE DISTANCE d(max 
AT WHICH THE CROSSWIND INTEGRATED CONCENTRATION AT GROUND LEVEI 
IS A MAXIMUM 


Data obtained at Idaho Falls, U.S.A., with a source at a height of 50 metres. 


Test Observed d (max) Calculated d (max) in metres 

300 370 

500 g10 

1000 850 

600 487 

500 445 

600 830 

11 400 305 
12 600 805, 
13 600 790 
14 400 540 
15 300 370 
16 700 765 
17 400 304 


The ratios of the calculated and observed values range from 0-76 to 1°82; 
apart from test 5 the range is 0-76 to 1°38; the overall average ratio is 1-14. 
These data suggest that useful estimates of d(max) may be made in a simple 
way from practicable measurements of wind fluctuations near the site of release. 


Vertical diffusion at longer range.— [he discussion of vertical diffusion 
has so far been restricted to distances of travel of about 1000 metres. For longer 
distances not only has there been no established treatment available, but until 
quite recently there were virtually no useful observational data. Consequently, 
estimates of diffusion at longer range have tended to be based on extrapolation 
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of the short-range data, a somewhat dubious procedure bearing in mind the 
primitive state of knowledge of turbulent and convective transport processes 
above the immediate surface layer. 


In the last few years tracer studies of diffusion over distances of tens of miles 
have been undertaken at Porton and elsewhere. The work at Porton has been 
partly concerned with horizontal spread, but the main interest has been to 
obtain some reliable description of the extent of vertical diffusion. Earlier stages 
of this work at Porton were dependent on a limited effort in sampling from 
aircraft, and this proved incapable of providing more than a very rough indi- 
cation of the vertical distribution of the tracer material. 


More recently a technique has been developed for sampling the cloud from 
a crosswind line-source of tracer material, using units mounted on a barrage- 
balloon cable, and this has given valuable preliminary data on vertical distri- 
bution at a distance of about 50 miles from a source. These data include a 
demonstration of relatively uniform concentration (with height) in a convective 
régime, with a relative sharp fall-off near the base of the overhead inversion (at 
about 3500 feet in the case studied). It was also found that in the absence of 
convection, but without any marked stabilization near the surface, vertical 
diffusion could be very slow. Two separate experiments showed the cloud to be 
essentially confined to the first 2000 feet above ground, and in one of these the 
material was actually released at about 1000 feet. This slow vertical spread is 
all the more noteworthy when it is realized that extrapolation from short-range 
data would lead to a cloud height of 11,000 feet. On the other hand, it is at any 
rate qualitatively consistent with the small vertical gustiness measured at the 
same time in the first few thousand feet. 


The latter measurements, with others made at intermediate distances, have 
been shown by Hay and Smith® to be consistent with a new statistical treatment 
of the spread of a cluster of particles (as distinct from a continuous plume). This 
treatment enables approximate estimates of the spread to be made, given merely 
the total intensity of turbulence (in effect the o, of equation (2) for large values 
of t), though for more detailed analysis a knowledge of the energy spectrum of 
the turbulence is required. 


A practical system for estimating the concentration or dosage pat- 
tern up to about 100 kilometres from a source (see Appendix I). The 
method set out in detail in Appendix I is an attempt to combine in the most 
flexible manner the various ideas and observations which are now available. 
The basic assumption is that the crosswind and vertical distribution in a plume 
or cloud of windborne material can be represented by the Gaussian form, as 
adopted in equation (3). This equation implies that the wind direction is steady 
over the duration of release or of sampling. For general application, including 
highly variable wind direction, it is more convenient to use a lateral distribution 
on an arc centred on the source (instead of a crosswind line). Assuming this arc 
distribution to be Gaussian, equation (7) in Appendix I follows directly from 
the continuity condition. 


Equation (7) gives the axial concentration from a ground-level source in 
terms of a lateral (angular) spread 6 and a vertical spread A, which are defined 
by concentrations one-tenth of the axial or ground values respectively. Remem- 
bering that for a Gaussian distribution these dimensions are respectively 4°3 
and 2-15 times the root-mean-square deviations from the axis or ground, it is 
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easily verified that equations (3) and (7) are identical when 4 is small, and z 
and H are set equa! to zero. 


The advantages of equation (7) are firstly its simplicity, and secondly its use 
of plume dimensions 9 and A which can be directly envisaged. There is accord- 
ingly less likelihood of unrealistic magnitudes of 6 and A being adopted. When 
the necessary special data on wind fluctuation are available, # should be cal- 
culated from equation (1). Likewise, it is recommended that 4 be calculated 
from equation (2), using data on the fluctuation of wind inclination well clear 
of the ground, except for short distances (say <1 km) from a ground-level source. 
For use in the latter circumstances, and also generally when data on wind 
fluctuation are not available, estimates of 4 and § in broad meteorological 
conditions are given in Figure 2 (Appendix I). 

The estimates of A in Figure 2 have the following origins: 

(i) Neutral conditions (D) and distance <1 km—experimental data con- 
solidated by Calder’s semi-theoretical treatment. 


(ii) Non-neutral conditions and distance <1 km—experimental data ob- 
tained in Project Prairie Grass in the U.S.A. For extremely unstable 
conditions a round figure of 1000 metres is adopted for A at a distance 
of one kilometre. 


For neutral—-moderately unstable conditions (D, C, B) and distance 
>1 km—calculations from the available statistics on vertical gustiness,® 
supported by some experimental data. 

(iv) Stable conditions (E and F) and distance > 1 km—these are essentially 
speculative extrapolations from the more reliable data. 


The estimates of @ tabulated on Figure 2 are for a short release (a few minutes 
and are based on recently acquired statistics of wind direction fluctuation."® 
Those for d=o-1 kilometre are derived from equation (1), while those for 
d=100 kilometres are extremely tentative values based on a little experimental 
data, and on Hay and Smith’s* treatment of an expanding cluster. 


Greater uncertainty in the data of Figure 2 is implied by the thinner and 
broken lines for A and the addition of brackets to the figures for 6. For longer 
releases no attempt is made to give statistical estimates of 6, and for use in the 
absence of detailed measurements of wind fluctuation a rough rule is given for 
deriving @ from a routine wind-direction trace (see Appendix I, para. 11). 

It will be noted that in the extension to an elevated source no attempt is 
made to adjust the values of vertical spread, A, (though this would be auto- 
matically introduced to some extent if A were calculated from equation (2) 
using measurements of wind inclination near the height of the source). This 
use of common values of A, irrespective of the height H of the source, may 
introduce additional error at short distances (or more strictly while h is small, 
and less than H say). At the present stage, however, there is insufficient data on 
which to base any general correction of practical consequence. It should also be 
emphasized that the estimates of 9 and A are appropriate to fairly level open 
country. In an urban area, or on an industrial site, there will be additional 
dynamical turbulence generated by the buildings, and this may be expected to 
increase the spread of the plume. No quantitative data are available, but a 
priori there would appear to be no reason to expect an important effect except at 
relatively short distances from the source, before the spread of the cloud becomes 
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large compared with the individual buildings. At these short distances buildings 
will also bring the further complication of downdraught and downwash effects, of 
the type described by Hawkins and Nonhebel," both of which tend to bring 
effluent to ground level more quickly than in unobstructed flow. 


The rest of the process of evaluating the distribution of concentration is dealt 
with in a self-contained fashion in Appendix I. It should be noted that no 
attempt is made to allow for deposition, decay or decomposition of the material. 


Acknowledgements.— The writer is indebted to Mr. N. E. Bowne and Dr. 
N. F. Islitzer, United States Weather Bureau, Idaho Falls, Idaho, for the data 
in Table I. 


REFERENCES 


. SUTTON, 0. G.; Micrometeorology. New York, London, Toronto, 1953. 

. HAY, J. 8. and PasguiLt, F.; Diffusion from a continuous source in relation to the spectrum 
and scale of turbulence. Advances in geophysics (Atmospheric diffusion and air pollution), 
New York, London, 6, 1959, p. 345. 

. CALDER, K. L.; Eddy diffusion and evaporation in flow over aerodynamically smooth and 
rough surfaces: a treatment based on laboratory laws of turblent flow with special 
reference to conditions in the lower atmosphere. Quart. 7. Mech., Oxford, 2, 1949, p. 153. 

. DEACON, E. L.; Vertical diffusion in the lowest layers of the atmosphere. Quart. 7. R. met. Soc., 
London, 75, 1949, p. 89. 

5. DEACON, E. L.; Wind profiles and shearing stress—an anomaly resolved. Quart. 7. R. met. Soc., 
London, 83, 1957, P- 537- 
MONIN, A. 8.; Smoke propagation in the surface layer of the atmosphere. Advances in geo- 
physics (Atmospheric diffusion and air pollution), New York, London, 6, 1959, p. 331 
7. CRAMER, H. E.; A practical method for estimating the dispersal of atmospheric contaminants. 
Proc. 1st nat. Conf. appl. Met., Hartford, Conn., Oct. 28-29, 1957, Boston, Mass., p. C-33. 

. SMITH, F. B. and Hay, J. s.; The expansion of clusters of particles in the atmosphere. Quart. 
J. R. met. Soc., London, 87, 1961, p. 82. 

. SMITH, F. B.; An analysis of vertical wind-fluctuations at heights between 500 and 5000 fect. 
(Awaiting publication. 

. SMITH, F. B. and ABBOTT, P. F.; Statistics of lateral gustiness at 16 metres above ground. 
(Awaiting publication). 

. HAWKINS, J. E. and NONHEBEL, G.; Chimneys and the dispersal of smoke. 7. Inst. Fuel, London, 


28, 1955, P- 53°- 


Appendix I 


Instructions for the estimation of the distribution of concentration 
or dosage downwind of a source of windborne material 


The Plume model and the formula for concentration or dosage from a ground-level source 


1. Consider a ground-level source producing a plume with an idealized 
distribution as represented in Figure 1. Let the lateral spread 6 ( =“ AOB) alongan 
arc be defined by concentrations one-tenth of the peak or axial value. Similarly, 
let the vertical spread A be defined by a concentration one-tenth of the ground 
value. With the simplifying assumptions that the wind speed u is constant with 
height, that the material crosses the arc normally and that the crosswind and 
vertical distributions are of Gaussian form, the concentration distribution is 
completely determined by equating the rate at which material crosses the arc 
to the rate of release at the source. 


2. For a rate of release of one “unit’’/min the axial or peak concentration C, 
at ground level at distance d downwind is 


2°8x10°3 


ud6h 


“units” /m4, 
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eh ad 
ONCENTRATION 


FIGURE I SCHEMATIC DIAGRAM OF A PLUME FROM A SOURCE AT GROUND LEVEI 


with u in m/sec, d in km, @ in degrees and A in metres. This is also the (otal 
dosage in “units” min/m* which would be experienced at the same position 
during the entire passage of the plume, when the (otal release is one unit. For 
any other rate of release, or total release, the values of C, should merely be 


increased in direct proportion. 


3. The formula is valid in terms of concentration and rate of release provided 
the duration of release is sufficient for alongwind diffusion to be neglected, and 
in practice this may be taken to be the case when the duration of release is, say, 
equal to or greater than the time of travel from the source O to the point ot 
interest P. Apart from a relatively short period near the beginning and the end, 
the concentration C, may then be regarded as obtaining (on average) for a 
period equal to the duration of release. On the other hand, when the release is 


terminated before the leading edge of the plume reaches the arc through P, the 
period of quasi-steady C, will be substantially less than the duration of release, 
and will be virtually nil when the duration of release is only a small fraction of 
the time of travel to P. In these cases, although equation (7) will give an over- 
estimate of the concentration experienced at P, it will still give the correct value 
of the total dosage for a total release of one unit. 


4. On either side of P, and vertically above any point on arc AB, the con- 
centration will fall off according to the (assumed) Gaussian form, and the 
complete downwind distribution may thus be determined from four para- 
meters: the speed of the wind u (appropriately defined), the appropriate 
effective wind direction (fixing the position P), the vertical spread h and the 
lateral spread 6. The wind values need to be based on surface and upper air 
data, in a way which is explained later. 


The estimation of the vertical spread h 


5. The magnitude of h initially increases with distance d from the source, at 
a rate which depends on the amount of vertical mixing. If, as frequently 
happens, vertical transport is suppressed at some level in the atmosphere by 
an isothermal or inversion layer, the ultimate effect of this should be to trans- 
form the concentration profile into a uniform distribution between ground and 
inversion (or isothermal) base. Beyond this stage the effective value of A will be 
constant. 


6. When data are available (either from current measurements or from 
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accumulated statistics) on the fluctuation of the wind inclination, values of A 
should be calculated from the appropriate form of equation (2), that is 

h = 2150d 0, metres, ..+ (8) 
where a, is the standard deviation of the wind inclination ¢ (radians), obtained 
from averages of ¢ taken cver periods equal to approximately one-quarter of 
the time of travel and observed for the duration of release or sampling.This 
procedure is recommended especially for an elevated source (for which con- 
dition corrections to equation (7) are made in para. 14), and also for a ground- 
level source once vertical spread has extended above the surface layer (10 metres 
or so in depth) in which wind shear is most pronounced. 


3000 


O03 O5 


morre 


FIGURE 2-——TENTATIVE ESTIMATES OF VERTICAL (h) AND LATERAL (9) sPREAD 


7. For use in the absence of wind fluctuation data, tentative estimates of 
vertical spread in open country are given in Figure 2 for six categories of 
stability (in the surface layer) which are specified qualitatively in terms of 
wind speed, insolation and state of sky. “Strong” insolation corresponds to 
sunny midday conditions in midsummer in England and “‘slight’’ insolation to 
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similar conditions in midwinter. Night refers to the period from one hour before 
sunset to one hour after dawn. The neutral category D should also be assumed, 
irrespective of wind speed, for overcast conditions during day or night, and for 
any sky conditions during the hour preceding or following night as defined 
above. The D(1) curve should be followed to the top of the dry adiabatic layer; 
thereafter, in sub-adiabatic conditions, D(2) or a curve parallel to D(2) should 


be followed. 


TABLE II-—KEY TO STABILITY CATEGORIES 
Insolation Night 
Surface ————, 


= 
wind speed Strong Moderate Slight Thinly overcast or <3/8 cloud 
at1om a » 4/8 low cloud 


EEE — 2 A. 





m/ $éC / 
A A-B B | 
A-B B Cc | E 
B B-C Cc D 
Cc Cc-D D D 
Cc D D D 








For A-B take average of figures for A and B, etc. 


8. In very light winds (<2 m/sec) on a clear night, that is, conditions pro- 
ductive of sharp ground frost or heavy dew, the vertical spread may be even 
less than the values given for category F. However, because of lack of quanti- 
tative knowledge of this and because in practice the surface plume is unlikely 
to have any definable travel, no estimates are attempted for this case. 


g. In unstable conditions the value of fh estimated as in paragraph 6 or 7 


/ 
should be used with increasing distance only until a magnitude A, is reached, 


equal to the estimated vertical limit of convection, at distance, say, d,. For 
approximate evaluation, at distances equal to or greater than about 2d,, a 
constant value of 2h, should be used. In calculating C, from equation (7) this 
allows roughly for the assumed development of the vertical distribution, from 
the Gaussian form at d, to uniformity at 2d,. (Exact allowance for a uniform 
vertical distribution over depth A, actually requires h=1-71 A, in equation(7). 


The estimation of the lateral spread 


10. When suitable data are available on the fluctuation of wind direction, 
the lateral spread should be calculated from equation (1). The 4 of equation 
(7) is equal to 4-3 times the oy of equation (1), when consistent units are used, 
the numerical factor being appropriate to the assumption of Gaussian distri- 
bution. 


11. When fine-structure data are not available, rough estimates of 6 for a 
long release, in the region of one hour or more, may be made from a routine 
wind-direction trace as follows: 


d = 0-1 km, difference between extreme maximum and minimum of trace 
over period of release, 
d = 100 km, difference between maximum and minimum “15-minute 
averages” of wind direction. 
For a short release (a few minutes) estimates are given on Figure 2 for the six 
stability categories. 





The evaluation of the axial concentration 

12. The axial concentration Cy for unit source strength may be calculated 
forthwith from equation (7) by substituting the distance d, the estimated values 
of h and 6 and the appropriate value of wind speed u. For most practical 
purposes it will be sufficient initially to do the calculation for four standard 
distances, 0-1, 1, 10, and 100 km. Moreover, since 4 usually changes only 
slowly with distance, it will be adequate to take values at 0-1 and 100 km, and 
to interpolate for 1 and 10 km by assuming equal changes in 9 in three intervals 
of distance. For the wind speed wu the ‘surface’ values should be used with 
d=o-1 and 1 km, but at the longer distances a mean value throughout the 
vertical extent of the plume is required. In practice, for vertical spread from a 
few hundred to say 1500 metres, a speed midway between the surface and geo- 
strophic speeds should be a reasonable working approximation. 


FIGURE 3—STANDARD VALUES OF AXIAL CONCENTRATION Cy FOR d=0-1 KM AND 
u==5 M/SEC 


For other values of d and u multiply by o- 1/d and 5/u (see equation (7 
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13. Rapid determination of C, may be carried out from Figure 3, which 
shows isopleths of C, for d=o-1 km and u=5 m/sec, for a practical range of h 
and 9. Corrections for distance and wind speed are easily made by multiplying 
the values from Figure 3 by o-1/d and by 5/u. The values for the four standard 
distances may then be plotted on log/log graph paper for subsequent inter- 
polation (Figure 4 shows such a graph for the example set out in Appendix 2). 
If the distance 2d, (see para. g) falls between 1 and 10 km, the line joining the 
points at 10 and 100 km should be produced backwards to 2d), and the inter- 
polated point so formed should be joined to the point at 1 km. If 2d, falls 
between 10 and 100 km, the line joining the points at 1 and 10 km should be 
produced forward to 2d,, and the interpolated point so formed joined to the 
point at 100 km. (in the example in Appendix 2, 2d, is 11 km but the difference 
between this and 10 km is disregarded, and the points at the standard distances 
are joined directly.) 


Allowance for elevation of source 

14. If the source is elevated the concentration at ground level will be reduced 
at the shorter distances but will tend more closely to the ground-source values 
as distance is increased, that is, as the cloud spreads vertically the initial effect 
of placing the source above the ground is progressively lost. Correction factors 
F, as a function of h/H,* where H is the height of the source, are given in the 
table below, and these should be used to multiply the previously derived 
values of Cy. For this purpose it will be necessary to take values of A at shorter 
intervals and this of course can be done using Figure 2. The quickest procedure 
is to evaluate the distances corresponding to the given values of h/H, and to 
apply the factors F, to the corresponding values of Cy on the graph of axial 
concentration (see Figure 4). 


A/H 1/2 2/3 4/5 1 1} 14 
F, 10-4 5°6x107% 0°027. O10) «©0°23~— 0 36 


15. In the special conditions referred to in paragraph 8 it is possible that at 
the height of an elevated source the wind speed may be sufficient to give 
appreciable travel of the plume, though vertical spread would be negligible. 
This would mean that the development of vertical spread should be started 
not at the source, but at a downwind position corresponding to the wind speed 
and the estimated time for breakdown of the stable situation. 


Plotting of the position and concentration of the plume (see Figure 5 for example 

16. In practice, interest will centre on the area covered by a concentration 
greater than some specified value. From the equivalent threshold value of 
concentration (that is, the actual threshold concentration divided by the 
source strength), the range of distance to be considered will be immediately 
evident. A drawing of the position and average distribution of the plume may 
then be prepared, as described below. If there is interest in the maximum 
distance of 100 km, then it will clearly not be practicable to construct the details 
within the 1 km distance, and for this a separate larger-scale drawing should 
be prepared. 





* From equation (3), F, = exp (-H*/2 oz? exp (-2°303H*/h?). 
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FIGURE 4——CONCENTRATION—DISTANCE DIAGRAM FOR THE EXAMPLE IN APPENDIX 2 





17. The first step is to estimate the axial positions P at the various distances 

and these are given by effective mean wind directions as follows: 

dinkm Effective wind direction 

Or 1-1 surface wind direction 

10-100 =average of the surface and geostrophic directions, backed by 10 
For a long release the basic estimates of wind direction should be made for 
the total period of release, and if direction changes with distance downwind 
appropriate allowance should be attempted from the synoptic data. Even in 
the case of a short release the important factor is the trajectory over some 
distance, and this will again be best given by an average wind direction cor- 
rected as necessary for variations downwind. Strictly speaking, the effective 
wind direction should be a mean through the vertical spread of the plume, 
weighted according to concentration. The rule given above for d=10 to 100 km 
is probably adequate for h of the order 1000 m, but for very much smaller or 
larger values increased weight should be given to the surface or upper winds. 


Ag 
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FIGURE 5-—PLAN OF PLUME FOR THE EXAMPLE IN APPENDIX 2 


18. Having marked the position O of the source, draw arcs at two or more 
standard radii (0-1, 1, 10 or 100 km). From O draw lines downwind, corre- 
sponding to the effective wind directions, to cut the appropriate arcs (at P, P, 
P, or P, respectively). The average plume axis is then given by joining OP, P, 
P, P, and values of concentration Cy can be entered on this as required, using 
figures interpolated from the graph of C, against d. On each arc mark off points 
A, B, etc. symmetrically about P, etc., so that the angle A, OB, is equal to the 
estimated lateral spread 9. Then OA, A, A, A, and OB, B, B, B, are the 
“boundaries” of the plume, at which the concentration falls to 1/10 of the 
corresponding value on the axis OP, P, P, P,. 


19. For positions other than on the axis or ““boundaries”’ concentrations may 
be interpolated by applying the following factors F,,* which are the ratios of 





* From equation (3), F, = exp/ »*/2e,? 





exp (2° 303(2a)?/6* 


47 





the axial concentration to the off-axis concentration, for deviation « from the 
axis. These factors enable the isopleth of a given concentration to be drawn in 
quickly, as follows. If the given concentration is C’, read off the concentration- 
distance graph (Figure 4) the distances at which the axial concentration is FC’, 
and plot F,C” on the axis of the distribution (Figure 5). Points on the isopleth 
of C’ are then obtained by marking off the corresponding deviations « from 
the axis. 
Deviation « from axis (fraction of 6/2 . as as tl! MUON 2 


) 
axial conc. 


. 18 2°8 4°4 10 37 180 104 


off-axis conc. 
Allowance for type of wind direction variation 
20. ‘The above procedure is based on a Gaussian distribution of the fluctu- 
ating wind direction, and this can be assumed to apply sufficiently well as long 
as the fluctuations are large compared with any discernible systematic trend 
over the period of release. On the other hand, if there is a systematic veering 
or backing of the direction, over a range large compared with the width of 
the trace, then the factors in paragraph 19 will not apply. The procedure 
should be exactly as before up to the construction of the plume position, but 
then the concentration should be taken as uniform along any given arc, and 
equal to 0-58C, (4C, will be adequate for most practical purposes), up to the 
cloud “boundaries” OA, A, A, A,. and OB, B, B, By. 


Rapid evluation of the distance and magnitude of the maximum ground-level concentration 
Srom an elevated source 

21. If interest or time is insufficient for carrying out the full procedure, an 
estimate may quickly be made of the distance d (max) along the axis, of the 
maximum ground-level concentration from an elevated source. Assuming that 
the lateral and vertical spread have the same variation with distance (that is, 
that 4/@d is independent of distance), d(max) is the distance at which Af is 
approximately 3/2 times the effective height H of the source. Having decided 
on the stability category, and given H, d(max) can thus be read directly from 
Figure 2. 


22. The magnitude of this maximum concentration, that is Cy(max), for a 
release of one “‘unit’’/min, is given by 
- Lo le 
C,(max) ~ ————_-__. “units’’/m*, paek@ 
3ud(max) 6H 
with u, d(max) and @ in the practical units adopted in equation (7), and H (the 
height of source) correspondingly in metres. 


Accuracy 

23. It is emphasized that the present system can in general give only very 
approximate estimates of the magnitudes of the concentrations, especially when 
it is necessary to use the tentative statistical estimates of A and 6. In the more 
difficult cases of unstable and stable situations, it is obvious that errors in A of 
several fold could be involved at the longer distances of travel, and this should 
be kept in mind in applying the data to the assessment of hazards. On the 
other hand there will be relatively straightforward cases when the estimates of 
vertical spread may be expected to be correct within a factor of two, namely: 
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PLATE U SUMMER SEA FOG SPREADI 
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AND BLOWING SNOW AT ADMIT 


IEPENDENCIES, MAY 1G 





PLATE IV HOLE IN STRATOCUMULUS, DECEPTION ISLAND, FALKLAND ISLANDS 
DEPENDENCIES, WINTER 1953 


The cloud was formed orographically by the ridges of Deception Islar 
} 


apparently produced by the freezing of a water-droplet clo and fil 


can be seen below the hole. There were no aircraft in the vicinit 





all stabilities except extremes, for distances of travel of a few hundred 
metres, in open country, 


neutral to moderately unstable conditions, for distances of a few kilo- 
metres, 


unstable conditions in the first 1000 metres above ground, with a 
marked inversion thereafter, for distances of travel of 10 kilometres or 
more. 


24. Uncertainties in the lateral spread of the plume are likely to be less 
important, except when the wind field is indefinite, in which case an even 
more important error will be that involved in prescribing the position of the 
plume. In such circumstances the best procedure would be to estimate the 
concentrations in the usual way, but then to allow for a wide range of possible 
directions of the plume, even to the extent of a full 360° in the most indefinite 
wind situations. 


Appendix 2 


Example of calculation of distribution of concentration from a point source of strength 
one “‘unit”’/min 


Figures in parentheses are relevant paragraph numbers in Appendix 1. 


Ge neral data 
Site 


Southern England 
Date 


Period 

Effective height of release (H 
Surface wind 

Geostrophic wind 

Vertical extent of convection 
State of sky 

Stability category (7 


16 June 1959 
1000-1300 GMT 
100 ™m 

4 m/sec, 275 

8 m/sec, 325 
1000 ™m 

1/8 Cu, 6/8 Sc, 
B-( 


Distance at which vertical spread (h), from Figure (2), equals ver- 
tical extent of convection : 5°5 km 


Calculation of Cy, equation 
Distance d . 100 km 
Effective value of A (7 2000 m 
Lateral spread 6 (11)* 

Effective value of u (12 ‘ d 6 


deg 
m/sec 
x J ° x . ‘3 units/m? 


Allowance for elevation of source (14 

Assumed h/H i/2 2 4/5 Oo 14 
F, exp 2°303H?/h? 10 5 « 3 9°027 0°10 0 

Distance at which correspond- 

ing values of A occur (Figure 2 0°28 0°37 
Gs * 4600 2700 1800 1200 580 3 
FC, 0°46 15 49 120 210 185 


2°0 4° 
36 0°56 0 


oO 46 0°59 0 86 1 F km 
2 
} 


© units/m?® 
© ynits/m?* 


Position of plume axis (18 

Distance hee id iad um ie 100 «km 

Effective wind direction (17 290 deg 

Rapid evaluation of d (max) and C, (max 
h at d (max 

“. d (max 

Ce (max 


21, 22 
—3H/2 150m 
from Figure 2 0°85 km 
from equation (g 2x 


3*~4 


107% 





xX 0°85 X95 x 100 


~ 2x10~ units/m® 
* Obtained from routine wind-direction trace. 
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DUST HAZE AT BAHRAIN 
By J. HOUSEMAN 


Poor visibility, caused by blowing dust and dust in suspension, is frequent over 
the deserts of Iraq and northern Arabia, especially during the summer months. 
The prevailing winds over these deserts are north-westerly. Bahrain Islands, in 
the Persian Gulf, lying to the south-east of the deserts, are consequently to 
leeward and are also affected by the dust. 


The increasing use of jet aircraft with their voracious appetites for fuel 
makes the forecasting of landing conditions more and more important. This is 
especially so for aircraft scheduled to land at Bahrain as, when dust is wide- 
spread over Iraq, Arabia and the Persian Gulf, the nearest unaffected airfield 
may be Damascus, over 800 miles away. An investigation into dust haze 
conditions was undertaken in the hope that it would be of value both to fore- 
casters and to aircraft operators for long-term flight-planning purposes. 


Reduction of visibility by dust at Bahrain Airport is partly influenced by the 
geographical position of the airfield and partly by local topography. The air- 
field is situated on the northern side of Muharraq Island, which lies north o/ 
Bahrain Island and some 25 to 30 miles east of the Arabian mainland. 


With only very few exceptions, dust-raising winds are of the north-westerly 
type, known locally as “Shamal”. Consequently dust affecting the airfield is 
nearly always brought from the Arabian coast, or from Iraq some 250 miles to 
the north-west. Locally lifted dust, unmixed with Arabian dust, while often raised 
over the main island to the south, only affects the airfield in squally conditions 
such as those connected with the passage of active fronts. It is infrequent over 
the airfield and is usually of short duration. 


In general, winds of at least 25 knots are needed to lift local dust at the 
airfield and winds of over 40 knots are required to give severe reductions in visi- 
bility. On the other hand, dust which has been lifted over Iraq or Arabia and is held 
in suspension can reduce visibility for days at a time with only very light local 
winds. This type of dust haze is therefore more frequent than that caused by 
dust raised locally and can be considered the main subject of the investigation, 
though it has not been possible to differentiate between the two forms as they 
occasionally occur in conjunction. 


Method.—The method used in the investigation is similar to that which has 
been used by a number of others in examining the frequency of fogs. The period 
investigated is from 1 January 1956 to 31 December 1958. The percentage 
frequencies of the occurrence of various visibilities for each hour of the day 
throughout the period have been evaluated and some of the results are illu- 
strated graphically in Figures 1 and 2. 


The three years show considerable variation in individual dust haze patterns 
and this variation obviously affects the results, but random sampling over ten 
other years indicates that the mean figures obtained are largely representative, 
provided only frequencies of 10 per cent and upwards are regarded as signi- 
ficant. The figures refer only to occasions when visibility was reduced solely 
by dust. Mist and fog are not included. 
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FIGURE I! PERCENTAGE FREQUENCIES OF VISIBILITIES 3 N. MILES 


Dust haze frequency and the climatological régime.— he main point 
of interest emerging from the figures is the way in which the dust frequency fits 
the climatological régime of the area. November to February is the rainy season 
and, although strong winds are then more frequent than in summer and the 
rainfall is slight, the moistening of the desert surface is still enough almost 
completely to prevent dust haze. In fact, in the period December to February 
dust was all caused by the passage of cold fronts accompanied by thunder- 
storms and squalls and was a mixture of locally raised dust and dust carried 
along by the fronts. 


During March the surlace soil begins to dry out and dust haze becomes more 
frequent, though the percentages are still low and the actual reduction in 
visibility is slight. Dust is partly brought from the mainland in suspension and 
is partly lifted locally by strong winds giving a small increase in frequency 
around midday when convection and gustiness are at their maximum. 


FIGURE 2——PERCENTAGE FREQUENCIES OF VISIBILITIES < 1100 YARDS 


From May, through June and July the dust is at its worst. This period, 
especially June and the beginning of July, is the time of the “‘forty-day Shamal”’, 
an almost continuous north-west wind of varying strength, usually about 15 to 
20 knots, associated with the monsoon low over India and Persia. It is notice- 
able that the highest frequency in these months occurs at 0400 Gat. This is 
caused by dust, picked up over Iraq during the afternoon and carried some 250 
to 300 miles at 15 to 20 knots, reaching Bahrain in the early part of the follow- 
ing day. 
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This tendency for the highest frequencies to occur in the early morning 
continues through August to October, although August, which is usually a 
month of light winds, often south-easterly, shows a very much reduced total 
frequency. In September the winter régime of short periods of strong north- 
westerlies alternating with short periods of lighter winds begins again, but with 
the first rains over the northern desert in October the dust-raising properties 
of the stronger winds are much curtailed and a rapid reduction in dust frequency 
occurs. The more general rain of November settles almost all the dust other than 
that raised by the occasional front. 


Occasions of very poor visibility.— Reduction of visibility below 2200 
yards is most infrequent and reductions below 1100 yards are even more so. 
During the period examined reductions below 1100 yards occurred only in the 
six months March, April, May, June, July and December and only in the 
month of June did they last for more than two hours. 

Average visibility.— Visibility in dust haze averages 3000 to 3500 yards. 
Once a slow improvement has brought visibility to three nautical miles, further 
improvement to six to ten nautical miles is usually rapid. During the winter 
months visibilities of 15 to 20 nautical miles are frequent. 

Conclusions.—From October to February the chances of aircraft move- 
ments being hampered by dust haze at Bahrain are almost negligible. March 
and April remain clear for night operations and are only occasionally hazy by 
day. During May, June and July dust may affect the airfield at any time, the 
worst month being June and the worst period being around dawn with a 
tendency towards improvement in the late afternoon. August and September 
are again usually clear during the afternoon and night but are occasionally 
hazy in the mornings, 


Even when dust is present average visibility is 3000 and 3500 yards. Reduc- 
tions to below 1100 yards are only likely in June or during the passage of a squall. 
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ROUTINE COMPUTATION OF MONTHLY UPPER AIR STATISTICS 
USING AN ELECTRONIC COMPUTER 
By D. DEWAR, B.Sc. 

Introduction.—In order to make full use of the results of the daily upper air 
soundings carried out for synoptic purposes at radio-sonde stations controlled by 
the Meteorological Office, it has been the policy of the Upper Air Section of the 
Climatological Division to have these data entered on suitable forms each 
month, and to use them for the production of monthly routine statistics giving 
values which are likely to be required at short notice for research, investiga- 
tions, or the supply of information for aviation requirements. 


In the early days of the Upper Air Section the computations were carried out 
by assistants with the aid of adding machines and tables only. Starting with 
data for 1948, much of the laborious arithmetic was eliminated by the use of 
punched cards and Hollerith machines but a lot of computing still had to be 
done by assistants using desk machines, and the checking of the computations 
and entries of final results on summary sheets required a considerable amount of 
assistants’ time each month. A further step forward has now been made, after 
some months of development work, by the use of machine “programmes” 
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coded instructions directing the operation of the computer) which enable 
nearly all the monthly routine computations to be done by the Ferranti 
Mercury computer installed in the Central Forecasting Office at Dunstable. 
A brief account of the procedure followed is given below. 


Checking of data and conversion to tape procedure.—<A _ rigorous 
checking system is used to try to ensure that errors in entries on the forms are 
eliminated before the forms are passed for the data to be punched on Hollerith 
cards and subsequently converted to the symbols on a reel of paper tape 


required for use in the computer. Errors which are not detected until the 
results of the machine processing are scrutinized entail an exorbitant waste of 
time in correcting the forms, Hollerith cards and data tapes. 


All the data required for the upper wind statistics for one ascent are given 
on one Hollerith card. Four Hollerith cards are used to record values of 
temperature, heights of isobaric surfaces and humidity at standard pressure 
levels for one ascent, and for technical reasons it was decided that it would be 
best to record on the tape the data from the first Hollerith card (giving data for 
the lowest four pressure levels) for each day of the month, then data for the 
second Hollerith card (giving data for the next four pressure levels) for each 
day and so on. Before conversion to tape, the cards are “‘sequence checked” on a 
collator; though not infallible this check virtually ensures that all the cards 
in a pack are those for the required station, year, month, etc. A further check 
is provided by hand-punching at the beginning of each data tape a “preamble”’ 
giving the station number, year, month and hour to which the data relate; 
instructions in the programme provide for these particulars to be checked 
against those given at the beginning of each daily card and if there is disagree- 
ment the computer prints REJECT followed by information which indicates 
why the data were rejected. ‘This checking procedure is considered worthwhile, 
although it adds to the operating time, as it also acts as a check on the correct 
functioning of the machines themselves. 


Development of machine programmes 


(a) Wind programme.—The development of the programmes was initiated by 
Mr. J. S. Sawyer, then Chief Forecasting Research Officer, who produced a 
demonstration “‘Autocode” programme for computing wind statistics using the 
daily wind components and also indicated the procedure to be followed to 
convert the data punched on Hollerith cards to a form, on tape, suitable for 
reading into the computer. Some modifications and additions have since been 
made to this programme to speed up the input of the data* and to provide for 
values for all hours combined to be given in addition to values for each hour of 
observation. Part of the print-out of results for Bahrain for October 1958 is 
shown in Table I. The programme provides for the output punch of the 
computer to give only essential indicative information and the required 
statistics; explanatory titles, numbering of columns, etc. have been added in 
italics in this illustration. In actual use, the results are filed in folders, one for 
each station, and a key to the data is provided inside each folder. 


Values of the “all hours” mean components, etc. are computed by dividing 
totals for all hours by the total number of observations as this practice had been 





* Mr. P. B. Sarson suggested this improvement and carried out the necessary amendments 
to the programme, 
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TABLE I-—-WINDS FOR BAHRAIN, OCTOBER 1958 


Combined hours data 
850 3°6 —0'5 5 3° ‘2 13 45 
700 3°8 ol 15°¢ 45 
500 . —13'0 290 3° . 18-: 44 
400 §°s 19°7 285 . 20° 43 
300 26-6 270 24° 4! 
200 —7° —29°9 256 a 26- 40 
150 =—12°6 —29°3 247 31° 35 26- 38 
100 -6-6 —12°0 241 3°7 6 15° 36 
a) Station number h) Vector mean wind speed 
6) Last two figures of year i) Scalar mean wind speed 
c) Month j) Standard vector deviation 
d) Hour of observation k) No. of observations 
¢) Mean N-S components 1) Direction of max. wind 
f) Mean E-W components m) Speed of max. wind 
(g) Vector mean wind direction n) Date of max. wind 
(p) Data for max. wind at any level—height (decametres), pressure, direction, speed, date 
Speeds are in knots; directions are in degrees from true north. 


followed in earlier years; a change to the method advocated by the Climato- 
logical Commission of the World Meteorological Organization in its Guide to 
climatological practices, whereby an “all hours” mean value is taken as the mean 
of the values for the separate hours, will probably be made in January 1961, the 
end of the current five-year climatological period. 


(6) Programme for temperature etc.—This programme is more involved than the 
wind programme as the data for one ascent are recorded on four Hollerith cards 
and three elements—temperature, heights of isobaric surfaces and humidity 
mixing ratios—are dealt with. 


Table II shows the print-out, for all levels for one hour of observation and 
for both hours combined, of statistics of heights of isobaric surfaces and 
temperatures. Column numbers in italics and a key to the values have been 
added but otherwise the values are as printed out from the output tape by a 
teleprinter. 


Statistics of temperature and height are first computed for the surface,* goo, 
850 and 800-millibar levels, and are then printed out and also stored for use 
in the “‘both hours” computations. This part of the programme is then repeated 
three times to deal with the other twelve pressure levels, four at a time. The 
programme contains instructions which cause the machine to print out days 
for which data are missing for a// the four levels being dealt with and to add 
in the number of these missing days to a count which it is instructed to make 
of consecutive days with missing data at any of the levels; if there are five 
or more consecutive missing observations, the machine prints an asterisk for 





* The height of the 1000-millibar surface is given in place of the surface height. 
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each such occurrence beneath the appropriate pressure-level figures. The 
object of this procedure is to allow a user to judge from the number of 
observations and the number of asterisks how much reliance he can place on 
values for that pressure level. At 50 millibars in Table II, for example, there 
were eleven observations and these, if there was no gap of five days or more, 
could just be regarded as giving a satisfactory mean value. The asterisks, 
however, indicate that there were two such gaps. It was not considered practi- 
cable to provide for this indication of reliability to be given for the “both hours” 
values also. 


The next part of the programme provides for the computation of similar 
values for both hours combined and for the differences between the 12 h and 
oo h means of temperature and height to be printed. These differences, 
in addition to providing useful information, enable a quick rough check of 
the mean values to be carried out. Provision has also been made for a rough 
check of extreme values to be carried out by the computer; if the range of 
either the isobaric heights or temperatures exceeds five times the appropriate 


standard deviation a query is printed in the space beneath the figures giving 
the difference between the mean values. 


TABLE II-—-TEMPERATURES AND HEIGHTS OF ISOBARIC SURFACES, FOR BAHRAIN, 


OCTOBER 1958 
a b ¢ 
40427 658 10 
e f g h j 
1000 2921 983 94 135 
31 31 21°4 63 
goo 31886 782 1029 1061 
31 31 17°4 ° 1003 
850 47339 668 1527 1555 
31 31 15°4 2 1505 
800 63481 2048 7 2079 
31 31 14°! 2029 


700 98236 3169 ; 3210 
31 31 13°8 3148 
600 137172 24 4425 ; 4466 
31 31 13°1 4402 
500 181860 242 5866 5.906 
31 31 19°3 5829 
400 234445 7563 - 7606 
31 31 26°7 7599 


300 298630 9633 9700 
31 31 36-2 , 9560 

250 337310 10881 10960 
31 43°0 | 10790 

200 6 12343 12430 
31 52°8 12230 

150 424120 14137 14220 
30 59°6 , 14030 


MISSING 16 
MISSING) 25 

100 479380 16530 - 16610 

29 49°5 *s 16450 

80 499050 17823 178g0 

28 39°5 17720 

60 410190 ; 19533 19620 

2! 44°7 19460 

227150 75 20650 20700 

11 37°9 . 20580 





TABLE IlI—TEMPERATURES AND HEIGHTS OF ISOBARIC SURFACES FOR BAHRAIN, 
OCTOBER 1958 (cont. 


é d 
10 BOTH 
f g 
5757 1836 
62 62 
63712 1583 
62 62 
946058 1395! 
62 62 
26088 1091 


b2 b2 


190550 505 
62 62 
274477 61 
62 62 
jo39911 472 
62 62 
jOogtio 1231 
62 62 


597010 215 


62 b2 


9560 
675000 2091 1o0ggo 
10790 
12470 
2220 


863030 39908 14270 


62 62 
705700 3309 


62 02 
oI o1 ‘ . 14000 


992620 1549 = 16543 3 = 166go 
bo 6o 50° 4 16430 
052370 $34! 17836 18000 
59 . 50 iF 4 17720 
)780g0 4 19558 é 19780 
50 ) 56-5 19460 
703560 ; 20682 20940 


34 7 69° 2 ; 20580 66 


l pper line Lower line 


Station number Height total No. of obs. 

Last two figures of year g) Temperature total No. of obs. 

Month Mean height S.D. of height 

Hour of observation Mean temperature S.D. of temperati 
e) Pressure level* j}) Max. height Min. height 

k) Max. temperature Min. temperature 
l) Difference 12h—ooh mean heights 
m) Difference 12h-ooh—mean temperatures 
Heights are in metres; temperatures are in “C. 

* Height values are for 1000 mb, temperatures for surface. 


rhe third part of the programme instructs the machine to compute daily 
values of humidity mixing ratio from the values of temperature and relative 
humidity (which were stored during input 


and values of saturation vapour 
pressure which were read in as part of the programme instructions. Statistics 
are then computed, printed and stored for the first hour and a similar procedure 
followed for data for the second hour. Data for both hours combined are then 
computed and printed. Results for Bahrain for. one hour and for both hours 
combined are shown in Table III. Explanatory headings have been added in 
italics but otherwise the Table is as printed out from the tape output. 


‘The “both hours’? mean values in this programme are obtained by taking 
the mean of the values for the separate hours. 
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III—HUMIDITY MIXING RATIOS FOR BAHRAIN, OCTOBER 1958 
H.M.R. DATA 


a 
40427 
e ; h 
1000 15 83 
goo 6°17 
850 2 4 6b 
800 3°69 
700 ] 2°57 
600 5 1°45 
500 20°¢ 0°07 
400 p 0° 28 


300 ’ 0°09 
$0427 58 BOTH 


1000 62 
goo b2 
b50 62 
400 8 be 
700 5Q° 62 
boo 62 
500 4 ; 62 
400 17° 62 
300 be 62 
Station number Pot 
Last two figures of year oO 
Month Mean 
Hour of observation Max. value 
Pressure level*® Min. value 
k) Difference 12h-ooh meat 
Humidity mixing ratio values are in gm ke~ 
* Values given against 1000 mb are actually surface 


Computation, output and printing-out times 

a) Wind programme Che initial reading-in of the programme instructions 
takes about one minute but this, of course, has only to be done once for all 
stations processed during a session. Reading in the data and computing 
statistics for one hour of observation takes about 23 seconds and punching out 
the results takes about 22 seconds. For a station making four ascents a day 
the total machine operating time is roughly 3} minutes. Subsequent printing 
out of the data from the output tape requires 9} minutes. 

(6) Programme for temperature etc.—The initial reading-in of this programme 
takes just over two minutes. Reading in and computing data for heights and 
temperatures for one set of levels (for example, 700, 600, 500, 400 millibars 


takes about 14 seconds and a further 17 seconds are required to punch out 


results. ‘1'o compute and punch out all the humidity mixing ratio values takes a 
little over a minute. The total operating time for a station making two ascents a 
day is a little over 6} minutes, but printing out the data from the output tape 
requires about 25 minutes. 


c) Total times.—Using the above test times and making allowances both for 
the time required to insert tapes in the tape reader and for shorter processing 
times for stations not making the normal two temperature and jour wind 
ascents a day, it is estimated that, with the present programmes, for the 24 or so 
stations to be dealt with each month, about four hours will be required to 
compute and punch out the data and about twelve hours to print out the 
values. It is hoped shortly to modify the programmes so as to be able to use 
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the line printer attached to the computer. This would greatly reduce the 
printing-out time and an output tape would then no longer be essential; if it is 
decided that this output tape could be dispensed with, it should be possible to 
reduce the total computing and printing-out time to something of the order of 
three hours. 

It is not possible, for several reasons, to give a satisfactory answer to the 
obvious question, how much time does the new procedure save ? Firstly, the 
work, to some extent, is still in the experimental stage and sufficient experience 
has not yet been gained to make a reliable estimate of times required for the 
ancillary processes—preparation of data tapes, scrutiny and, if necessary, 
correction of results necessitated by incorrect basic values or faulty printing. 
Secondly, the computer does far more than could be done by the assistant staff 
available. One interesting comparison can however be given. Some years ago 
similar wind statistics to those obtained from the computer were worked up by 
assistants from Hollerith tabulations, though only for ‘all hours combined’. 
The production of the required Hollerith tabulations took more time than is 
now required for the conversion of the cards to tape for use by the computer. 
The computation and checking of statistics worked from the Hollerith tabula- 
tions for the 20 or so stations then dealt with took a time equivalent to that of 
one assistant’s work for about 15 days. Using the computer, a reasonable time 
per station for computing and punching out results for “both hours combined” 
is about 130 seconds and for printing out results a little under two minutes is 
required, that is, about 14 hours for 24 stations. 


The best answer probably is that the new procedure permits many more 
data to be computed using fewer staff for this side of the work and, what is 
perhaps more important still, it substitutes intelligent scrutiny by the assistants 


for the repetitive numerical drudgery they were formerly required to carry out. 


REVIEWS 


Atlantic hurricanes, by G. E. Dunn and B.-I. Miller gin. x 6in., pp. xx +326, 
illus., Interscience Publishers Inc., 250 Fifth Avenue, New York 1, 1960. 
Price: $10. 


This is a very good book. It is probably the most comprehensive book which 
has yet been written on hurricanes and covers every aspect of the subject (except 
mathematical theories). It should certainly be read by every forecaster in the 
area and by anyone contemplating research on the subject. There are eye- 
witness accounts from the ground and from the air, details of the life-history of 
hurricanes and a comprehensive summary of techniques for forecasting move- 
ment and development though, as the authors agree, there is still much room 
for research in this field. A chapter, which would be of particular interest to 
readers living in areas affected by hurricanes, concerns preparations to guard 
against damage. 


The book is well printed, and the tables and diagrams are well placed in 
relation to the text. Most of the diagrams are clear, though the numerals in 
those on pages 176 and 177 are too small. A friend of the reviewer, not a 
professional meteorologist, describes it as ‘‘a most readable book’’ and adds “‘a 
word of praise must be given to the very useful Glossary of Meteorological 
Terms and the Indexes”. 


5. E. VIRGO 





Magyarorszdc, Eghajlati Atlasza, Klimaatlas Von Ungarn (Climatological atlas of 
Hungary). 19 in. x 13 in., pp. 20+ 78, illus., Akadémiai Kiad6-Budapest, 1960. 


The atlas under review contains 130 charts most of which are on a scale of 
1:1,250,000 and the rest on 1:2,500,000. The contents and descriptions are in 
German as well as Hungarian on a loose inset. 


The volume opens with four maps of topography, soil types and vegetation. 
The climatic charts cover hours of sunshine, cloudiness, fog, actual (that is, 
unreduced) temperature, frost duration, dates of first and last frost, frequencies 
of “‘frost’’, “ice”, “extreme heat’’ etc. days, vapour pressure, relative humidity 
at 1400, precipitation amounts and frequencies of days of precipitation of over 
one millimetre etc., evapotranspiration, amount of snowfall, duration of days of 
snow and snow lying, first and last dates of snow, wind roses and mean isobars, 
synoptic charts for six characteristic weather situations, temperature extremes 
shown by mean isotherms of the warmest and coldest winter and summer 
months, extremes of precipitation in similar form and phenology. Pillar 
diagrams of mean monthly rainfall at a number of stations are also given. The 
periods vary but are mostly for 1901 to 1950 and two of the extreme months 
are taken from before 1910. The charts for the various elements are not all 
given for the same portions of the year; the portions are changed to suit the 
natural variations of the elements. No element has twelve-monthly mean charts 
but means for mid-season months are given for many. The selection has been 
carefully made to give the best impression of the variations and extremes in 
minimum space. 


The printing is very good. The areas between isopleths are coloured to 
harmonious schemes which vary with the element. This atlas is an excellent 
production and the reviewer’s only criticism is that he would have liked to 
have seen German as well as Hungarian titles and legends on the maps 
themselves. 


G. A. BULI 


Weather forecasting for aeronautics, by J. J. George. 10} in. x6} in., pp. ix+673, 
illus., Academic Press, New York and London, 1960. Price: £5 7s. 6d. 


The title of this book gives a good idea of its contents except that the words 
“in the U.S.A.” should have been added. The author, assisted by seven other 
contributors, is a meteorologist with Eastern Air Lines, and the book is almost 
wholly concerned with weather forecasting in relation to the operation of an 
airline in the United States of America. 


In the first somewhat philosophical chapter the author stresses the economic 
value of weather forecasts. No forecast is certain; each forecast has a probability 
of success, and when considering the economic operations of an airline it is the 
probability of occurrence of a weather phenomenon which must be evaluated. 
The simple economics equation is that protective measures for an eventuality 
should be taken if the probability of its occurrence P > C/L where C is the cost 
of taking preventative measures, and L is the loss which would result if these 
measures were not taken and the adverse eventuality materialized. Prediction 
diagrams are extensively used in the forecasting techniques described in the book 
and from these, in most cases, the probability of the meteorological eventuality 
can be estimated. 
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Eight chapters are concerned with the production of forecast charts. The 
chapter on the prediction of cyclogenesis (45 pages) is very thorough: the author 
describes methods of determining whether or not cyclogenesis will occur, and if 
it does whether it will take the form of deepening of the parent cyclone, or 
formation of a new cyclone, or of a centre jump. He then describes objective 
methods of determining the location of new cyclogenesis, its timing, the future 
intensity, the future track and the predicted speed. An even more extensive 
chapter (94 pages) is concerned with the movement, deepening and filling of 
cyclones, and objective methods are described for forecasting their development, 
speed, direction of movement and time of recurving. 


These two are the most comprehensive chapters, but very useful and impor- 
tant are the following five chapters dealing with the movement of anticyclones 
in North America, the movement of cold lows at the 500-millibar level, the 
displacement of surface cold fronts, warm frontal analysis and movement, and 
the movement of tropical cyclones. A chapter entitled “the poor-man’s numerical 
weather prediction system” completes the chapters concerned with the produc- 
tion of forecast charts. 


The next four chapters are concerned with forecasting weather phenomena. 
There is a chapter on the prediction of very low ceilings and fogs, followed by 
one on pre-trough winter precipitation in which the author produces prediction 
diagrams for forecasting the amounts of precipitation. There then follows a 
chapter on the prediction of severe weather, thunderstorms, line squalls, tur- 
bulence, hail, tornadoes and aircraft icing—the latter is dismissed in a page as 
a phenomenon which can no longer be regarded as a hazard. Heavy snowstorms 
are discussed in a short chapter, wind and temperature forecasting at length, 
and there is a very lucid short chapter on the use of radar. The last hundred 
pages of the book are concerned with local forecast studies for airfields in the 
United States. 


This is a first-class stimulating book. It is generously illustrated and it is so 
well written that the formidable task of reading 662 pages becomes at once a 
pleasure. What is wholly admirable about the book is its straightforwardness: 
it is what it professes to be, a handbook for forecasters, and while underlying 
theories are sometimes mentioned the emphasis is in describing forecasting 
methods. What is impressive is the variety of the forecasting methods; the 
parameters for the prediction diagrams are derived variously from surface, 
850, 700 and 500-millibar charts, sometimes they are isotherm gradients, some- 
times wavelengths or amplitudes; clearly the methods result from a great deal 
of trial and error, based on theory, and the charts and methods described are 
those which have been found by experience to give the best results. 


The main value to the British forecaster of reading the book lies in its stimulus 
since probably none of the forecasting methods described as suitable for the 
United States can be used in the United Kingdom without modification. The 
book shows what can be done if a real effort is made. A great deal of operational 
research has been carried out in the United States of the kind that is necessary 
to bridge the gap between the theoretical work of the Rossbys, the Scherhags 
and the Sutcliffes, and the operational demands of the forecaster on the bench. 
This excellent book might be regarded as part of the dividends paid by that 
operational research. 


R. A. HAMILTON 





OBITUARY 


Mr. Edward William Barlow, B.Sc.—It is with deep regret that we record the 
death of Mr. E. W. Barlow in his 74th year, on g January 1961, after a retire- 
ment of only two years. A short sketch of his career appeared in the Meteoro- 
logical Magazine of January 1959. 


Mr. Barlow’s 39 years’ service in the Meteorological Office, which he joined 
as a Senior Professional Assistant after service with the R.N.A.S. and R.A.F. in 
the First World War, was noteworthy in that he spent no fewer than 32 years 
in the Marine Division. Here he was primarily concerned with ocean currents 
and sea-ice, the preparation of the atlases on these subjects, based on observa- 
tions sent in by British ships, and the corresponding text in the 73 volumes of the 
Admiralty Pilots, all of which he revised once, a number twice and a few three 
times, as fresh knowledge made a new edition possible. 


He was responsible also for the selection of items from the meteorological 
logbooks of ships of the voluntary observing fleet for publication in our con- 
temporary, Zhe Marine Observer, and as an expert on various natural pheno- 
mena frequently contributed explanatory notes for them. His particular study 
was bio-luminescence, and on this subject he continued to advise the Marine 
Division after his retirement. L.B.P. 


HONOUR 
The following award was announced in the New Year Honours List, 1961: 
C.B.E. 


W. A. Grinsted, Director of the West Indies Meteorological Service. 


METEOROLOGICAL OFFICE NEWS 


Retirements.— The Director-General records his appreciation of the services 
of: 


Mr. C. 7. G. Budd, Senior Experimental Officer, who retired on 28 December 
1960. He joined the Office in June 1920 as a Technical Assistant at Croydon. 
In 1925 he was transferred to the Aviation Services Division at Headquarters 
where he remained for some nine years, except for a short spell at an aviation 
outstation in 1927. Since 1934 he served continuously at aviation outstations 
including a tour of duty at Malta. From 1948 until his retirement he served at 
Uxbridge. 


Mr. H. Gingell, Senior Assistant (Scientific), who retired on 17 December 
1960. He joined the Office in January 1936 as an Observer, Grade II, The 
greater part of his service has been spent at aviation outstations including a tour 
of duty in Iraq. He also served for short spells in the Instruments Division in 
1945-6 and the British Climatology Division in 1954. In 1955 he was transferred 
to the London Forecasting Office where he remained until his retirement. 


Staff suggestions scheme 

Mr. F. B. Swain, Experimental Officer, was awarded {£25 for a suggestion 
leading to the introduction in the Meteorological Office of simplified topo- 
graphical maps. 
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TELECOMMUNICATIONS 


speed weather forecasting 


Plessey is everywhere, you'll find ... in the Middle 
and Far East for example, where minute by minute plotting 
and data transmission are the only key to effective warning 
of adverse weather conditions. 
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The ubiquity of Plessey Radio Teleprinter Equipment is as complete in 
world meteorology, press and consular services, as Plessey Telecommuni- 
cations is in the total telecommunications field. A vital amenity for any 
civilised community, the accurate dissemination of weather data is 
essential to airways, to planned agriculture, and in many localities to the 
very preservation of life itself. It is a field of service in which the Plessey 
Company is proud to be making a major contribution. 

Behind this achievement are the extensive prototype and quality 
manufacturing resources of the Telecommunications Division, inspired by 
the work of a unique research and development organisation. The 
Company recognises the need to segregate advanced thinking from the 
hustle and distraction of the factory. Accordingly, in addition to the 
Group Research Laboratories at Roke Manor and other specialised re- 
search centres already existing in the U.K., extremely well equipped 
laboratories have recently been established at West Leigh for advanced 
telecommunications studies. In these establishments, the next generation 
of telecommunications equipment is already taking shape. 


In close support at all times are the complete resources of the Plessey 
Group of Companies which include unrivalled tool making and machining 
facilities, a full range of environmental testing and production laboratories, 
and the services of the materials laboratories at Caswell. 


TELECOMMUNICATIONS DIVISION - ELECTRONIC & EQUIPMENT GROUP 
THE PLESSEY COMPANY LIMITED: ILFORD - ESSEX ~- Telephone: ILFORD 3040 
Overseas Sales Organisation: Plessey international Limited « Ilford * Essex * England 
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Observatory 


Founded 1877 


A Magazine presenting current developments 
in Astronomy by means of Articles, Corres- 


pondence, Notes on Discoveries and Reviews 


of important astronomical books. The papers 


read at the Meetings of the Royal Astronomi- 
cal Society and the discussions which follow are 


also fully reported. 


Annual Subscription for 6 issues, 


40s. 
including postage 


should be sent to 


“The Observatory’? Magazine 
Royal Greenwich Observatory 
Herstmonceux Castle 
Hailsham, Sussex 
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the most economical windfinding system 


An outstanding feature of the new Decca WF2 Windfinding Radar is its high 
performance coupled with low daily operating cost. A new high power transmitter 
increases the possible maximum slant range to some 200,000 metres, and enables a 
smaller and cheaper target to be used for shorter ranges. The equipment is so designed 
that it can be operated by one man and is capable of providing wind speed and direction 
up to an altitude of 100,000 fect in all parts of the world. 


WF2—the most advanced windfinding radar in the worid 


iAH Be CIEE OECCA RADAR LIMITED - LONDON - ENGLAND 











SELECTED PUBLICATIONS 
OF THE METEOROLOGICAL OFFICE 


A Century of London Weather. By W. A. L. Marshall (M.O. 508, 1952.) 


15s. {post 9d.) 
Condensation Trails from Aircraft. (M.O. 479, 2nd edition, 1952, reprinted 1954.) 
Od (post 2d.) 


Handbook of Aviation Meteorology. (M.O. 630, 1960.) (post Is. 9d.) 


Handbook of Meteorological Instruments: Part 1. Instruments fos 


surface obser- 
vations. (M.O. 577, 1956.) 


455 (post ls. 9d.) 


Handbook of Weather Messages, Codes and Specifications. (M.O. 510.) 
Part I, Transmission schedules and station index numbers. (3rd edition, 1959.) 
8s. (post 7d.) 
8s. (post 7d.) 

Part Ill, Coding, decoding and plotting. (3rd edition, 1959.) Bs 
(Amendments issued as necessar 


Part Il. Codes and specifications. (3rd edition, 1959.) 
(post 7d.) 


y and priced separately) 


Obtainable from 
HER MAJESTY’S STATIONERY OFFICE 


or through any bookseller 











TECHNICAL OFFICERS 
required by the 


EAST AFRICA HIGH COMMISSION 


for the Meteorological Department, mainly for weather forecasting, 
on contract for one tour of 30/36 months in first instance. 


Commenc*ng salary (including Inducement Pay) according to 
experience in scale £1,056 rising to £1,956 a year. Gratuity at rate 
134 per cent. of total salary drawn. Outfit allowance £30. Free 
passages. Liberal leave on full salary. 

Candidates must have Higher School Certificate with Physics and 
Mathematics or equivalent qualifications. Successful 


candidates 
without forecasting experience will be required to undergo course 


of training lasting 21 weeks before leaving U.K. in July. 


Apply to CROWN AGENTS, 4 Millbank, London, S.W.1 for 
application form and further particulars, stating age, name, brief 


details of qualifications and experience and quoting reference 
M3C/53164/MI 
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